Quantum phase transitions in d-wave superconductors 
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Motivated by the strong, low temperature damping of nodal quasiparticles observed in some 
cuprate superconductors, we study quantum phase transitions in d^2_y2 superconductors with a 
spin-singlet, zero momentum, fermion bilinear order parameter. We present a complete, group- 
theoretic classification of such transitions into 7 distinct cases (including cases with nematic order) 
and analyze fiuctuations by the renormalization group. We find that only 2, the transitions to 
d^2_y2 -\- is and d^2_y2 + id^y pairing, possess stable fixed points with universal damping of nodal 
quasiparticles; the latter leaves the gapped quasiparticles along (1,0), (0, 1) essentially undamped. 
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Recent photoemission [ij and Thz conductivity Q mea- 
surements on Bi2Sr2CaCu2 08+5, the cuprate supercon- 
ductor, have indicated anomalously large inelastic scat- 
tering of fermionic quasiparticles near the gap nodes 
in the d-wave superconductor. While many scattering 
mechanisms and scenarios have been proposed 0, 0) B SI 
for the damping of quasiparticles along the (1,0), (0, 1) 
directions (the "anti-nodal quasiparticles" ) above the su- 
perconducting critical temperature Tc, the possibilities 
below Tc at the nodal points are much more restricted, 
and allow us to make sharp distinctions between compet- 
ing theories. Standard BCS theory predicts a nodal scat- 
tering rate ~ from short-range interactions, and this 
is far too small to account for the observations. In this 
paper we study a possible explanation^ due to proximity 
to a quantum phase transition to some other supercon- 
ducting state X (see Fig [IJ . We show how global sym- 
metry and field-theoretic considerations permit a classi- 
fication of all possibilities for X, and we list those that 
may account for the experiments. 

The nodal quasiparticles at the gap nodes have a mo- 
mentum distribution curve (MDC) with a width propor- 
tional to ksT and there is little changeQ] in this be- 
havior when tuning T through Tc. The anti-nodal quasi- 
particles are broad and ill-defined above Tc, but narrow 
dramatically below Tc, forming long-lived states with an 
energy gap of 30-40 meV. A natural possibility, based on 
other experimental probes Q , is that these effects are due 
to proximity to a quantum critical point to magnetic or- 
dering. However, wavevector matching conditions appear 
to rule this out for the nodal quasiparticles: the magnetic 
fluctuations are strongest near wavevector Q = (7r,7r), 
and while they can strongly scatter anti-nodal quasipar- 
ticles above Tc, they do not connect low energy quasipar- 
ticles near the nodes 

Rather than exploring the intricate details of the many 
experiments, this paper performs the following well- 
posed theoretical task: classify and describe theories in 
which a d-wave superconductor at[l^ T <^Tc has, with 
minimal fine-tuning, (a) a nodal quasiparticle MDC with 
a width oc fc^T, and possibly {b) negligible scattering of 
the quasiparticles along (1,0), (0,1). We find that the- 
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FIG. 1: Phase diagram adapted from Reflrl. Superconduc- 
tivity is present for T < Tc. The long-range order associated 
with the state X vanishes for T > Tx , but fiuctuations of this 
order provide anomalous damping of the nodal quasiparticles 
in the quantum-critical region. 



ories which satisfy (a) also have a high frequency tail 
ll( in the energy distribution curve (EDC) of the nodal 
quasiparticles, as is experimentally observed P, 0|. 

Strong scattering of the gapless nodal quasiparticles 
surely requires their coupling to some low-energy bosonic 
mode. It is convenient to imagine that we have at our 
disposal some parameter r (which is possibly the hole 
concentration 6, but not necessarily so) which we can 
tune to condense the bosonic mode, leading to a new su- 
perconducting state X for r < rc (Fig[T|). The quantum- 
critical region of the phase transition at r = rc and T = 
will satisfy (a) provided the phase transition is below its 
upper critical dimension, and the nodal fermions arc in- 
trinsic (in a sense to be made precise below) degrees of 
freedom of the critical field theory [l^]- Conversely, (6) 
requires that the anti-nodal fermions are merely specta- 
tors of the phase transitions, and are essentially decou- 
pled from the critical degrees of freedom. 

The most efRcient scattering of nodal quasiparticles 
is provided by a linear, non-derivative coupling between 
the fermion bilinears and the order parameter; higher 
order and derivative couplings have been considered re- 
cently 0, [3], and invariably lead 0] to quasiparticle 
scattering rates that vanish with super-linear powers of 



2 





E d 


2C4 21 21' 


basis functions 


s 


1 1 


1 


1 1 


1 


p 


2 -2 








(singa:, singy) 




1 1 


-1 


1 -1 


cos qx — cos qy 


dxy 


1 1 


-1 


-1 1 


sin qx sin qy 


9 


1 1 


1 


-1 -1 


sin qx sin qy (cos qx — cos qy) 



FIG. 2: Character table of the irreducible representations of 
the group dv The C4 rotations are about the z axis, and 
the / (/') are reflections about the (1, 0), (0, 1) ((1, 1), (1, -1)) 
directions; the basis functions are chosen to be invariant under 
translations by reciprocal lattice vectors. 



T. Order parameters which carry a net momentum Q, 
will, by momentum conservation, couple linearly with 
the nodal fermions only if the spacing between two of 
the nodal points is exactly Q. Transitions involving the 
onset of spin 14 1 or site/bond charge density waves 
(stripes) do satisfy 01 {a,b) for such values of Q; how- 
ever the restriction on Q could be a fine-tuning condi- 
tion, and is not satisfied by the Q values observed so 
far. "Staggered-flux" orderp, jl3l| has a derivative cou- 
pling to the nodal fermions, and Q — (tt, tt) which does 
not connect nodal points: so (a) is not satisfied. In- 
deed, only the value Q = naturally satisfies the con- 
straints of momentum conservation, and so we limit our 
attention to order parameters at zero momentum. Fur- 
thermore, spin-triplet ordering at Q = implies ferro- 
magnetic correlations which are unlikely to be present, 
and therefore we further restrict to spin-singlet fermion 
bilinears. This means that our order parameter is a com- 
ponent of the complex superconducting pairing function 



(cq|C_q|), or the real excitonic (or 'particle-hole') 



pairing function 



^qa'-^qa/ 



(cqa annihilates an elec- 



tron with momentum q and spin a It is useful 

to decompose the functions Aq and Aq into components 
which transform under one of the irreducible representa- 
tions of the symmetry group of the Hamiltonian 15| : this 
is Ciy X Z2 , where 6*4^ is the tetragonal point group (see 
Fig [2]), and the Z2 component represents time-reversal 
symmetry T (point group symmetry breaking has been 
considered recently [l| 



17[ J as have exciton condensations 
[13| | at non-zero Q). Generically, a second-order transi- 
tion can only occur by condensation of an irreducible 
component (multiple components can appear in succes- 
sive transitions), and this permits a complete classifica- 
tion of inequivalent order parameters. Note that dx2_y2 
pairing is already present for r > (see Fig [T]) , and 
we will assume that this ordering remains well-formed 
across the transition; all our subsequent characterizations 
of possible orderings in state X refer to additional order- 
ing beyond an implicitly assumed background of d^2_y2 
pairing. is necessarily even under T, and so can gen- 
erate s,p, . . . exciton ordering; similarly Aq can generate 
s,p, . . . pairing or is, ip . . . pairing (the latter also break 



T), leading to a total of 15 possible order parameters 
for X. Of these, s exciton ordering is equivalent to an 
innocuous shift in the chemical potential, while p and 
ip pairing are forbidden by Fermi statistics. Because of 
the background dx2-y2 pairing, further dx2_y2 or idx2-y2 
pairing is not a new ordering, while simple symmetry 
considerations {e.g. examination of the fermion disper- 
sion relation in state X) show that g excitons, g pairing, 
and d^2_y2 excitons are equivalent to dxy pairing, dxy 
excitons, and s pairing respectively. Only 7 inequivalent 
order parameters now remain and we will discuss their 
properties shortly. 

We begin by reviewing the action for low energy 
fermionic excitations in a d-wave superconductor. We 
denote the components of Cq^ in the vicinity of the four 
nodal points (iif, ±K) [K « 0.397r at optimal dop- 
ing) by (anti-clockwise) /la, /2a, ha, ha, and intro- 
duce the 4-component Nambu spinors ^E"! — {fia, ^abfli,) 

and *2 = {ha,£abflb) where £ab = Sba and e^i = 1. 
Expanding to linear order in gradients from the nodal 
points, we obtain 

r (pi^ , 

= -7^-^T'^^\{~iuJn + VFkxT'' +VAkyT'')^i 

/(P k 
-— 2 *2 (-*^« + VpkyT' + VAkxT"") ^2- (1) 

Here u;„ is a Matsubara frequency, r" are Pauli matrices 
which act in the fermionic particle-hole space, kx^y mea- 
sure the wavevector from the nodal points and have been 
rotated by 45 degrees from qx.y co-ordinates in Fig[21 and 
Ui?, f A are velocities. 

We now describe the 7 possible order parameters for 
state X, along with the respective actions for the quan- 
tum phase transition. 

(A) is pairing: This has been considered in Ref 0- 
The state X (with dx-^-y^ + is pairing) has no gapless 
fermionic excitations, breaks T, but all charge neutral 
observables (like the charge density or lattice displace- 
ments) retain the full C41, symmetry. The order parame- 
ter transforms as a real, one-dimensional representation 
of Ciy X Z2, and so can be represented by a single, real 
field (j)] this will also be true for (B)-(F) below, with only 
(G) requiring a doublet of real fields. On general sym- 
metry grounds, following action for 4> is obtained after 
integrating out high energy fermion modes: 



d^xdr 



i(a.0)^ + y(V0)^ + ^'^^ + ^'^'' 



(2) 



here r is imaginary time, c is a velocity, r tunes the sys- 
tem across the quantum critical point, and u is a quartic 
self-interaction. By itself, would describe a critical 
point at r = Tc in the universality class of the classical, 
three-dimensional Ising model. However, a coupling to 
the low energy fermionic modes in ([1]) can preempt this 
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conclusion 1J|: its form can be deduced from the val- 
ues of the basis function in Fig [2] at the nodal points, 
and the information that the order parameter is in the 
particle-particle channel- 



while the coupling between 



and Vti 2 is 



d?xdT 



(3) 



where A is the required linear coupling constant be- 
tween the order parameter and a fermion bilinear, and 
Ml = Ma = tV. 

(B) idxy pairing: This is very similar to A, with the main 
change arising from the new basis function in Fig[21 which 
now implies Mi = —M2 = . 

(C) ig pairing: Also related to (A), but now the basis 
function in Fig [5] vanishes at the nodal points. Conse- 
quently, the coupling between 4" and </> requires at least 
one spatial derivative, and is irrelevant The action 

in ([2]) is the entire critical theory of the transition, 
and the scattering of the nodal fermions is weak, arising 
only from irrelevant couplings, and violates (a). 

(D) s pairing: T remains unbroken, but the symmetry 
of charge neutral observables is broken to C2V , so that X 
(with dx'^-y^ + s pairing) is a superconducting nematic 
[3, [l3|- The nematic order is polarized along the (1,0) 
or (0, 1) directions. For weak ordering, the state X re- 
tains gapless nodal fermionic excitations, but the nodal 
points are at (±if', ±K) with K' ^ K; for a sufficiently 
large s component, the nodal points disappear upon col- 
liding in pairs as min(_ftr', K) — > 0, in a separate quantum 
critical point which is not of interest here. As in (A,B), 
coupling of the order parameter is described by Q, but 
with Ml = M2 = . 

(E) dxy excitons: This is as in (D), but symmetry of 
charge neutral observables in X is broken to a different 
C2V subgroup of Ciiv, with the nematic now polarized 
along the diagonal (1,±1) directions. The nodal points 
in X are at ±{K,K) and ±{K\ -K') with K ^ K'. In 
the action ([3]), we now have Mi — — M2 = r^. 

(F) dxy pairing: Such an ordering in X moves the nodal 
points clockwise (or anti-clockwise) from {±K,±K), 
reducing the C'^y symmetry to C4, while preserving 
T. Again the action ([3]) describes the order parame- 
ter/fermion coupling, but with Mi = — M2 = t^. 

(G) p excitons: The order parameter transforms under a 
two-dimensional representation of C^y, requiring a dou- 
blet of real fields, {(j)x,(j)y), to describe the low energy 
bosonic modes. The state X retains T and the gapless 
nodal points, but has broken to Z2. The action ^ 
is replaced by 



Sfj, = J d^xdr 

+ 4 {dy(t>.f 



Ci{dx4>yf + C\{dy(t>yf + e{dx(t>x){dy(t>y) 



d^xdr 



(5) 



(4) 



We now make a few general remarks on the field the- 
ories above. Upon integrating out the fermion fields, we 
find a finite one-loop renormalization of the tuning pa- 
rameter, r. This should be contrasted with the behavior 
in a system with a Fermi surface, where we would find 
the BCS infrared logarithmic divergence in the analo- 
gous term: this is, of course, the reason that a T = 
Fermi liquid is unstable to superconductivity for any at- 
tractive interaction. In the present situation, the back- 
ground dx2_y2 superconductivity has reduced the Fermi 
surface to 4 Fermi points, and so further pairing or ex- 
citonic instabilities occur at finite values of r and A. In- 
deed, this feature allows a non-trivial quantum critical 
point, with a universal quantum-critical region (Fig[T]); 
the fluctuations in this region will satisfy (a) provided the 
quantum-critical point at r = r^, T = is described by 
a fixed point of the renormalization group (RG) trans- 
formation at which A approaches a non-zero and finite 
fixed point value — then the scattering rate of the nodal 
fermions will be determined by T alone [isj . 

The results of our RG analysis of (A-G) are simple 
and remarkable. Only for (A,B,C) do we find a fixed 
point, accessed by tuning the parameter r; such a fixed 
point describes a second-order quantum phase transition 
at the critical point r = Vc- For all other cases, we 
find runaway flows of the couplings, with no non-trivial 
fixed points, which suggests first-order transitions. As we 
have already noted, the fixed point for (C) is the Ising 
model — the nodal fermions are decoupled from the crit- 
ical degrees of freedom in the scaling limit, so that (a) 
is not satisfied. Only (A) and (B) satisfy (a), with the 
couplings A and u approaching non-zero fixed point val- 
ues: the nodal fermions and (j) are strongly coupled in 
the critical theory, and the anomalous dimension of the 
fermion field leads to a large w tail in its EDO 0, El]. 
The (A, B) fixed points are also Lorentz invariant — the 
dynamic exponent z — 1, and the velocities renormal- 
ize to vp — = c in the scaling limit. Indeed, these 
fixed points were discussed earlier [7| , but only for almost 
equal velocities; here we have established that the equal- 
velocity fixed point is the only one for arbitrary initial 
velocities. However, the crossover exponent which deter- 
mines how rapidly the velocities approach each other is 
negligible (« 0.05), so that a transient regime with 
unequal velocities will be realized over essentially all of 
the experimentally accessible regime. 

The methodology of our RG is standard and details 
appear elsewhere 19]. The familiar momentum-shell 
method, in which degrees of freedom with momenta be- 
tween A and A — liA are successively integrated out, fails 
here: the combination of momentum dependent renor- 
malizations at one loop, the direction-dependent veloci- 
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ties {vf, f a, c . . . ), and the hard momentum cutoff gener- 
ate unphysical non-analytic terms in the effective action. 
So we obtained the RG equations by using a soft cutoff 
at scale A, and by taking a A derivative of the renormal- 
ized vertices and self energies. We obtained equations 
for all the velocities, the dynamic exponent z, and the 
field anomalous dimensions to one-loop order in the non- 
linearities A, u, V. For (D,E,F) a simple and robust effect 
preempts a fixed point: the structure of Mi ^2 produces 
opposite sign renormalizations for w_f,a, in a manner that 
both flow equations cannot simultaneously be at a fixed 
point; (G) required a more detailed analysis. 

Our main result is that, among the 7 transitions con- 
sidered here, only for those involving onset of d^^-y^ +is 
or d^2_y2 + idxy pairing in a dxi_y2 superconductor did 
we find a universal critical theory of coupled fermionic 
and bosonic order parameter modes below its upper crit- 
ical dimension. Such transitions naturally satisfy (a). 
Upon further imposing condition (&), case (B), with 
dx2-y2 -\-idxy pairing, is uniquely selected: from the basis 
functions in Fig [T] we see that couples to fermions in all 
directions for (A), while the fermionic coupling vanishes 
along the anti-nodal directions for (B) — so the gapped 
anti-nodal fermions will [will not] lose the sharp quasi- 
particle peak below Tc by emission of multiple quanta, 
for (A) [(B)]. 

Pairing in the d^^-y^ + id^y channel has been consid- 
ered in numerous works recently (20j . with the order in 
the ground state either global (induced spontaneously or 
by an external magnetic field) or local (in the vicinity of 
defects 21 1, surfaces [2^ . or vortices 23]). Here we only 
require strong fluctuations of such order, induced by a 
proximity to a hypothetical point in the phase diagram 
where global order arises. While experimental discov- 
ery of such a point is of course preferable, tests of our 
proposal would also be provided by signals of <j) fluctu- 
ations. This is a spin-singlet mode with dxy symmetry, 
odd under time-reversal, and at T = it has spectral 
weight with mean frequency and width both of order of 
an energy scale ~ {r — rcY'^ (where v is the usual correla- 
tion length exponent) — we estimate this scale is ~ 5 — 10 
K; in the quantum-critical region the characteristic en- 
ergy scale is kBT/h. Fluctuations of should lead to 
anomalies in Raman scattering [24] and Hall transport 
[2^: these issues will be discussed in future work. 
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Erratum: Quantum phase transitions in (i-wave superconductors 
Phys. Rev. Lett. 85, 4940 (2000) 

Matthias Vojta 

Institut fiir Theoretische Physik, Universitdt zu Koln, Ziilpicher Str. 77, 50937 Koln, Germany 

Ying Zhang 

Goldman Sachs, 1 New York Plaza, New York, NY 10004 
Subir Sachdev 

Department of Physics, Harvard University, Gamhridge, MA 02138 

We correct an error in our paper Phys. Rev. Lett. 85, 4940 (2000) [arXiv:cond-mat/0007170| . 
Our characterization of the physical properties of the superconducting state G was incorrect: it 
breaks time-reversal symmetry, carries spontaneous currents, and possesses Fermi surface pockets. 



In the discussion of case G in Ref.[li above Eq. (4), the 
single sentence "The state X retains T and the gapless 
nodal points, but has C^^ broken to Z2" is incorrect. The 
state X — G breaks T (time- reversal), and has sponta- 
neous electrical currents. For (j)^ and (j)y = (or vice 
versa) the currents have the same symmetry as those in 
the state On discussed by Simon and Varma 0. Also, 
as pointed out by Berg et al. the nodal quasiparti- 
cles do not survive in the superconducting state G, but 
turn into Fermi pockets. The latter conclusion can be 
verified from the fermion spectrum obtained by diago- 
nalizing Eqs. (l)-|-(5) for constant (px.y- 

All other sentences and the conclusions in the paper 
[H remain unchanged. 

Also, in the companion paper, Ref. 4, the only error is 
in the sketch of the fermion excitations in Fig. 2 for case 



G. 
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